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Single polymer adsorption in shear: flattening versus hydrody- 
namic lift and surface potential corrugation effects 
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PACS 82.35.Gh - Polymers on surfaces; adhesion 
PACS 83 . 50 . -v - Deformation and flow 

PACS 82 . 37 . Gk - STM and AFM manipulations of a single molecule 

Abstract. - The adsorption of a single polymer to a flat surface in shear is investigated using 
Brownian hydrodynamics simulations and scaling arguments. Competing effects are disentangled: 
in the absence of hydrodynamic interactions, shear drag flattens the chain and thus enhances 
adsorption. Hydrodynamic lift on the other hand gives rise to long-ranged repulsion from the 
surface which preempts the surface-adsorbed state via a discontinuous desorption transition, in 
agreement with theoretical arguments. Chain flattening is dominated by hydrodynamic lift, so 
overall, shear flow weakens the adsorption of flexible polymers. Surface friction due to small- 
wavelength surface potential corrugations is argued to weaken the surface attraction as well. 
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The adsorption of polymers on surfaces is at the base 
of many applications such as surface modification, col- 
loidal stabilization and flocculation. By tailoring adsorp- 
ion properties, polymer additives can be used as lubri- 
cants, adhesives, strength enhancing agents or coagulants. 
Many experimental and theoretical investigations on the 
basic research level exist, see [TJ for references. The the- 
oretical foundation for polymer adsorption has been laid 
down by de Gennes [2]. Recent works fine-tuned the mi- 
croscopic picture and addressed the adsorption of charged 
or stiff polymers [3j[4] . 

Although it was realized early that non-equilibrium ef- 
fects are very relevant for polymer adsorption because of 
the long relaxation times |2j, most theoretical and exper- 
imental work concentrated on equilibrium aspects [l][5]- 
However, technological applications involving polymer ad- 
sorption are typically far from equilibrium. Along the 
same lines, many biological mechanisms involve adsorp- 
tion of macromolecules in shear flow, e.g. the initiation 
of the coagulation cascade involving the van-Willebrand- 
factor [6| or the adsorption of E. coli on surfaces [7]. Non- 
equilibrium aspects of polymer adsorption receive growing 
attention both from the experimental [8 15 and theoreti- 
cal point of view 16 -22 but are still less well understood 



than the equilibrium case. In the context of the current 
theoretical investigation, two opposing effects are relevant: 
A polymer that is either pulled laterally by a terminally 



ing surface is flattened. In the absence of hydrodynamic 
interactions, this has been found to enhance adsorption 
in simulations 
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On the other hand, hydrodynamic 
effects generate a shear-dependent long-ranged lift force 
that decays quadratically with distance from the surface 
for dumbbell models 23 as well as for stiff and flexible 
polymers [24JI26] which in simulations has been found to 
weaken polymer adsorption 20 . In the present work we 



exerted force 25 or subject to shear 15 22 at an adsorb- 



disentangle those opposing aspects and consider flexible 
polymer at adsorbing surfaces in shear flow using simula- 
tions and theory: In the absence of hydrodynamics, we ra- 
tionalize the shear-induced flattening of the chain and the 
adsorption enhancement observed in simulations by linear- 
response theory. As we turn on hydrodynamic interactions 
in simulations, long-ranged hydrodynamic lift effects dom- 
inate over chain flattening effects and lead to chain des- 
orption via a discontinuous transitions, as predicted by 
scaling arguments. An analysis of a single-particle model 
based on the Fokker-Planck equation suggests that lateral 
motion over a surface with an inhomogeneous or corru- 
gated potential weakens the effective attraction, in line 
with previous simulations |21| Our overall conclusion is 
that shear flow weakens the adsorption of flexible polymer 
chains on planar surfaces. 

In the simulations, a polymer is composed of N spherical 
beads with diameter a. Neglecting particle inertia, the 
time evolution of the position of monomer i obeys the 
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Fig. 1: a) Schematic picture of the polymer model for non-equilibrium adsorption in shear. The z axis coincides with the 
surface normal. The position of monomer i is denoted as Vi = {xi,yi,Zi), the distance between two monomers i and j as 
Tij — (xij ,yij , Zij) — rj — ri. The force F m is the force exerted by monomer m + 1 on monomer m. A simple shear flow in 
x direction acts with shear rate y. b) The mean height profile {zi/a}, i = 1 . . . 128 as a function of the monomer index i for 
a 128mer obtained from equilibrium BD simulations at the adsorption strength x = 2 and decay length k _1 = 1/2. c) The 
mean rescaled and averaged lateral force Fi/j — {F^^fio/ja as a function of the bond index i for a 128mer. The profiles are 
obtained from the height distributions shown in b) via eq. Q. 
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Ti(t + At) -fi(t) = -^2 frijV? .«({f at}) 

+ ^2 Vfj-fiij + ikilhex + &(*)• (4) 



a typical setup of the simulation is shown in fig. [T|a). 
Hydrodynamic effects are incorporated via the mobility 
matrix fiij which is obtained from the Green's function 
G w (rj,rj) of the Stokes equation that satisfies the no- 
slip condition on a planar boundary |26||28] . Performing a 
multipole expansion to second order in terms of the bead 
diameter a, which is accurate for bead-bead and bead- 
surface separations larger than w 2a, we write 29 
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for i ^ j. The self mobility tensor fj, u is obtained in 
the limit — > Vj and regularized far away from the 
surface and set equal to the bulk sphere mobility fiQ = 
1 / {Z'Krja) |29| . This hydrodynamic treatment is valid 
in the incompressible low-Reynolds-number limit and as- 
sumes instantaneous propagation of hydrodynamic effects, 
which is accurate for most experimentally relevant situa- 
tions. Effects due to the discreteness of water molecules 
are in fact negligible for all but sub-nanometer length 
scales |30j. In the free-draining simulations, the mobility 
tensor is set to be diagonal and constant in space, 



(3) 



where 1 is the diagonal unitary tensor. The externally 
imposed flow is a linear shear with shear rate 7. For the 
Brownian dynamics simulations, we use the rescaled, dis- 



All lengths are scaled by the bead diameter, r = r/a 
and energies by thermal energy, u({rjv})fc,B J 1 = C/({rjv}) 
which leads to rescaled mobilities p,ij = n^AthsT/a 2 
with timestep At and to a rescaled shear rate 7 = 
7a 2 /(ksTfio). The random velocities which couple the 
system to a heat bath arc modeled with Gaussian white 
noise and fulfill the fluctuation dissipation theorem, 
(£i(t)£j(t')) = 2ftij5(t — If). The rescaled time step 
fid = fioAtkBT/a 2 must be chosen small enough such that 
the bead displacement per time step is small compared to 
the bead radius. The total potential u = u sur f + u po i y 
consists of bead-bead interactions 
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where 77, = \fj — fi\ is the rescaled monomer distance. 
The first term ensures the chain connectivity by har- 
monic bonds around the equilibrium length a with a 
rescaled spring constant k — 200, the second is a trun- 
cated Lennard- Jones potential with a rescaled parameter 
e = 2 which is only used for fij < 1 to avoid overlap of the 
chain. For the attraction between surface and monomers 
an exponentially decaying potential is used, 



^surf 
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where 



L /a and x are the rescaled decay length 



and the rescaled interaction parameter, respectively. For 
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Fig. 2: Free draining (FD) results, a) Simulation results for the mean normalized adsorption potential — («surf}/A r x for a 
128mer as a function of the surface interaction parameter \ at various shear rates 7. b) Results from the linear-response theory 
in which the lateral stretching force distribution due to shear is calculated from the equilibrium monomer height distribution via 
eq. |7| and used to derive the effective adsorption strength via eq. ( 10 I, see text for details. The linear response theory captures 
the shear-induced adsorption enhancement only qualitatively, c) Fraction of adsorbed monomers (T) /TV from simulations. A 
monomer i is defined adsorbed if Zi < ft . 



electrostatically driven adsorption, the interaction param- 
eter x can be interpreted as the product of surface and 
polymer charge density and the decay length kT 1 as the 
screening length. We use the exponential potential as a 
generic form to study short-ranged surface attraction and 
have in mind also hydrophobic or van-der-Waals interac- 
tions. To make it very short-ranged, we use k _1 = 1/2. 
We also include a hard wall interaction of the monomers 
at the wall-liquid boundary, i.e. u SUI {(zi < 0) = 00. 

We use a rescaled time step of /2o = 10~ 4 and simu- 
late for at least 10 s simulation steps. The first 20% of 
these steps are disregarded for equilibration. Every 10 3 
steps the configuration is recorded. Data are obtained by 
block-averaging, statistical errors are obtained from the 
standard deviation of blocks and only shown when larger 
than the symbol size. Simulations arc performed for poly- 
mers of typical length N = 64 or 128 and for different 
values of the adsorption strength (x = 10~ 4 . . . 10 2 ). 

We first briefly review the adsorption of a single poly- 
mer in equilibrium. In thermodynamic equilibrium, any 
polymer of finite length will desorb from an adsorbing 
surface into the semi-infinite half-space, no matter how 
strong the adsorption potential is. However, by carefully 
choosing polymer length and simulation duration, infor- 
mation about polymer adsorption in the asymptotic limit 
of an infinitely long polymer can be obtained. This is so 
because a finite-time window exists within which an ad- 
sorbed polymer of finite length has enough time to equili- 
brate its conformation but at the same time stays trapped 
in the surface adsorption potential. In fact, this time win- 
dow widens with increasing monomer number N, or, con- 
versely, shrinks for too short polymers. 

The transition between the adsorbed and desorbed state 
can be quantified with different observables, e.g. the ad- 
sorption potential energy [5j, the polymer mean height 



[4] or the number of adsorbed monomers 16 . In fig 
[2^ a) we plot the mean adsorption potential per monomer, 
— ( u smf/xN), calculated from eq. ([6|, for a 128mer as a 



function of \. In the thermodynamic limit N — > 00 this 
observables is expected to go to zero in a continuous fash- 
ion at the adsorption transition. 

Next we establish a simple linear response theory for 
the flattening of a chain in shear flow. For simplicity, 
we consider a freely jointed chain (FJC) and neglect hy- 
drodynamic interactions and excluded-volume interaction 
between beads. We assume the chain to take some average 
separation from the surface. From eq. (TlJ) the x-component 

(x) (x) 

of the velocity of bead i reads Xi — jzi+noF,- — //o-F/lj + 

(x) (x) 

Q where F- denotes the force that bead i + 1 exerts on 

(x) 

bead i via the bonding potential. By reciprocity, —F^_[ is 
the force that bead i — 1 exerts on bead i and the bound- 



ary condition reads as F 



(x) _ jp(x) _ 



JV 



0. By averaging 



this equation, the random force disappears and we obtain 
(ii) = j(zi)+no(F^) — ^q(F^\). In the stationary state, 
all monomers have the same average velocity equal to the 
mean chain velocity, {x\) — (± 2 ) = ■ ■ ■ , which furnishes 
N — 1 equations. Solving for the forces, we obtain 
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Mo 



(7) 



fc=i 



which constitutes an exact relation between the average 
chain height profile and the mean forces acting on the 
polymer bonds. 

According to standard polymer statistical mechanics, 
the mean extension of a FJC bond along the pulling direc- 
tion is (xu+i) = (xi+i— Xi) — C{Fj) where C(x) = cothx— 



x~ is the Langevin function and Fi = a{F> )/ksT is a 
short-hand notation for the rescaled average x-component 
of the bond force. Likewise , th e second moment is given 



by 



(i 2 ) 



1 - 2C(F l )/F l 
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If we neglect effects of 



the surface, the two perpendicular directions, y and z, are 
equivalent and the mean squared extension perpendicular 
to the pulling direction is 



(il) = {l-(ii))/2 = C{F i )/F i 



(8) 
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Fig. 3: Mean normalized adsorption potential —{u SUI {) /N\ 
for a 64mer as a function of the surface interaction parame- 
ter x in equilibrium (solid lines) and for various shear rates 
7 = 7a 2 / fioksT. a) Results from free-draining (FD) BD sim- 
ulations and b) from hydrodynamic (HI) simulations. In con- 
trast to the FD the HI simulations show weaker adsorption and 
a discontinuous desorption transition. 



with the asymptotic limits (zf) F -i^r 1/3 for weak force 

and (zf) 1/1 -Ft I for strong force. Since the root mean 
squared bond radius defines the Kuhn length, we can in 
the presence of the lateral pulling force define an effective 
Kuhn length in z direction by a^ ff ; = 3a 2 £(Fi)/ Fj. In 
fact, since the force varies along the chain contour, it is 
useful to define the renormalized Kuhn length as an aver- 
age over all bonds, 



N - 



Fj 



(9) 



For illustration we show in fig. [IJb-c) the height profile 
and the force profile at x = 2. Due to the symmetry of the 
mean height profile (zi) of the chain and since the chain 
ends have a larger separation from the surface than the 
chain middle, according to eq. (m), the force vanishes on 
the middle and terminal bonds and its magnitude is max- 
imal in between. Note that a non-zero bond force profile 
would even arise for e.g. polymer loops where by symme- 
try the mean height profile is fiat, caused by chain fluctu- 
ation around the mean (a mechanism not pursued further 
in this paper). Based on the mapping of the problem of a 
single polymer at an adsorbing surface onto the equivalent 



problem of a quantum particle at an adsorbing wall 32 



out to be the effective adsorption strength 



(10) 



This gives a direct clue as to how the chain flattening 
enhances chain adsorption. 

As exemplary case we show in fig. [2]ja) the mean nor- 
malized adsorption potential — {u sm f)/(Nx) for a 128mer 
from free-draining simulations as a function of the adsorp- 
tion strength x for the equilibrium case (solid line) and for 
different shear rates 7 (symbols) . For low values of the ad- 
sorption strength \ we note shear-enhanced adsorption in 
agreement with previous free-draining simulations 17p9 



At large values of \ the adsorption behaviour is not mod- 
ified. Intuitively, one might expect peeling and tumbling 
effects to reduce adsorption at strong shearing, but noth- 
ing of this sort is seen. In fig. [2]; we show the fraction of 
adsorbed monomers or trains, (T) /N, where a monomer 
counts as adsorbed when its height is less than the adsorp- 



tion screening length, 



< K 



in the strong adsorption 



regime. Almost no influence of shear is seen. This is eas- 
ily understood when taking the effects of chain flattening 
into account. For high adsorption strengths, the poly- 
mer adopts a flat configuration with long trains and only 
few, small loops. The stretching forces acting on the poly- 
mer bonds are low and chain flattening is negligible and 
thus the enhanced adsorption mechanism discussed previ- 
ously does not come into play. On the other hand, close 
to the desorption transition loops proliferate and experi- 
ence stronger shear forces, which leads to chain flattening 
and thus enhanced polymer adsorption. On the linear- 
response level, the configuration and adsorption energy 
of a chain in shear can be approximated by the equilib- 
rium configuration of a chain (i.e. without shear) with the 



the only remaining scaling variable in the problem turns 



adsorption strength x replaced by XeS given by eq. (10). 
The result of this procedure is shown in fig [2jb) . The 
linear prediction qualitatively captures the shear-induced 
adsorption enhancement but grossly overestimates the ef- 
fect. In fact, linear response breaks down at strong shear 
flows since the resulting chain flattening of the chain very 
efficiently decreases the shear effects, this is even more en- 
hanced by the stronger adsorption of the flattened chain. 
The linear theory can in principle be self-consistently im- 
proved, which we do not pursure since in the following we 
show that hydrodynamic lift effects dominate chain flat- 
tening. 

In fig. [3] we plot the rescaled adsorption potential 
— ("surf )/^Vx as a function of the adsorption strength x 
for the equilibrium case and various shear rates for chains 
of length N — 64. In a) we show the results from free- 
draining (FD) simulations, in b) we show data for the ex- 
act same parameters using full hydrodynamic interactions 
(HI). As expected, the equilibrium results (solid lines) 
from both methods coincide within the precision of the 
simulations. For moderate shear rates the HI simulations 
quantitatively confirm the results obtained in the FD case. 
For elevated shear rates, however, the transition into the 
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Fig. 4: (a) Potential energy u of a single particle at fixed height z = in a surface potential eq. (Ill and pulled with 
force of F = —0.1. The surface interaction strength is chosen as y = 1, corrugation power is p = 6. (b) Average 



/i//io for a single particle pulled through a corrugated potential eq. (Ill as a function of the rescaled pulling force F - 
for fixed particle height I = 0, /? = 6, k = 1, x = 6.28 and x = 62.8. Data points are taken from simulations, the 
Fokker-Planck theory eq. \V2\ . (c) Average surface potential — (ii su rf) for a single particle as a function of its surface 
z from stationary solutions of the Fokker-Planck eq. (121 with surface potential parameters h = 1, j3 = 6, and x = 
intermediate force F = 100 the average potential crosses over from the equilibrium limit [F = 0, solid line) to the 
limit (F = 1000, dashed line) at intermediate separation. 



a rescaled 
mobilities 
= Fa/k B T 
lines from 
separation 
: 62.8. For 
high-force 



desorbed state changes from continuous (FD) to discon- 
tinuous (HI) and the magnitude of the adsorption energy 
is smaller compared to the equilibrium case. We conclude 
that hydrodynamic effects weaken adsorption in shear and 
totally dominate the chain flattening effects seen in the 
absence of correct hydrodynamic effects. 

What is the cause for this drastic hydrodynamic repul- 
sion? It is known that dumbbells 
rods and flexible polymers 
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as well as stiff 
close to a flat, non- 
adsorbing surface in shear experience a repulsive lift force 
away from the surface. This lift force can be rationalzed by 
the anisotropy of the hydrodynamic mobility parallel and 
normal to an elongated object [27] in conjunction with 
the anisotropically distributed orientation of the object 
in shear [26]. The resulting potential is long-ranged and 
decays as oc 1/z for rods as well as flexible polymers at 
elevated shear rates [26] . The combination of such a long- 
ranged repulsive potential with a short-ranged attractive 
surface adsorption potential is known to turn the adsorp- 
tion transition of a polymer discontinuous 



33 . Our sim- 



ulations confirm the predicted change in the nature of the 
transition. In summary, due to the relative weak enhance- 
ment of adsorption due to chain flattening observed in 
free-draining simulations, the hydrodynamic lift repulsion 
dominates and weakens the adsorption and actually turns 
the transition discontinuous. 

Up to now the adsorbing wall was assumed flat and 
homogeneous. This is a good approximation for atomi- 
cally flat surfaces or adsorption driven by hydrophobic in- 
teractions since here the surface friction effects are small 
Also for charged substrates the assumption of a 
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laterally homogeneous potential is accurate if the screen- 
ing length is larger than the monomer bond length and 
the separation between charged sites on the substrate. In 
all other cases, and in specific for adsorption driven by 
hydrogen-bonds 35 , adsorbing polymers do experience 



corrugated potentials the non-equilbrium effects of which 
will be briefly discussed now. Since the equilibration of a 
polymer in a corrugated potential landscape including hy- 
drodynamic interactions is by far too demanding from the 
computer time aspect, we consider the simplistic case of 
a single monomer dragged by an external force in a one- 
dimensional corrugated potential. As a straightforward 
extension of the potential eq. ^ used in the simulations, 
we consider the total potential u(x, z) = u sur f (x, z) — Fx 
including the corrugated surface potential 



u sur{ (x,z) = -xexp(-«5) 



cos (2irx) + 1 



(11) 



where the exponent /3 controls the corrugation steepness, 
the corrugation wavelength is for simplicity set equal to 
the monomer diameter a, and the force F leads to par- 
ticle motion. Here we concentrate on the highly corru- 
gated case P = 6 shown in fig. Qa) for z = 0, F = —0.1 
and x = !• We consider numerical solutions of the sta- 
tionary Fokker-Planck (FP) equation 36 = P(x, t) = 
u'(x)P'(x,i) + u"{x)P{x,i) + P"(x,t), where P{x,t) is 
the normalized distribution function, the rescaled time is 
t = t^oksT/a 2 and u'(x) = du(x)/dx denotes a partial 
derivative w.r.t. x (and similarly for P). Note that z only 
enters parametrically into the FP equation, i.e., the parti- 
cle is for illustrative purposes constrained to fixed height 
z and we therefore dropped the z-dependence of the po- 
tential. For time-independent potentials consisting of a 
linear and a part periodic over the interval from x = to 
x = 1, the stationary solution is given by 



36 



P(x) 



-u(x) 



J,, 1 dx'e u{S: "> 
1 _ e «(i)-«(o) 



dx'e u{i ' ] 



(12) 



where Q is the normalization constant. The 
rescaled mobility follows as fi/fio = (x)/(F£to) = 
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- J dxu'(x)P{x)/{FQ) and the average potential energy 

follows as (u sur f) = Jq 1 dxu sm { (x)P(x)/Q. In fig. 
[4jb) the normalized particle mobility fi/fio from the 
FP solution (shown as lines) compares not surprisingly 
favorably with the Brownian dynamics simulation results 
(data points) according to eqs. ^ and Q at fixed 
particle height z — 0. Data is presented for \ — 6.28 
(circles) and \ — 62.8 (crosses). For low external force, 
the particle mostly sticks to the attractive surface sites 
and the average mobility is very small. At large forces 
the average mobility approaches its bulk value and the 
motion is uniform regardless of the potential corrugation. 
In fig. |4|c) we show the average potential energy — (u SU rf ) 
as a function of the vertical height z for fixed surface 
interaction parameter \ = 62.8 and for different driving 
forces. For comparison, the equilibrium case (i.e. for 
vanishing force) is shown as a solid line. For the larger 
force F — 1000 the mobility is not modified (as seen from 
fig. |4^b) for z = 0) by the potential and the particle 
moves uniformly over the surface, the average potential 
is therefore a uniform average over the corrugation. In 
contrast, in equilibrium, the particle mostly samples 
the potential minima and the magnitude of the mean 
potential is much larger. For intermediate force F = 100 
an interesting phenomenon appears: For small values 
of z the particle sticks to the potential minima and the 
average potential is close to the equilibrium case. With 
increasing height z the lateral confinement due to the 
potential is weakened, the mobility increases, and the 
average potential crosses over to the high-force limit. 
This means that the effective potential in the presence of 
a lateral driving force becomes more short-ranged, under 
the influence of shear where the driving force increases 
with distance from the surface this effects will be even 
more pronounced. 

For polymer adsorption in shear this means that po- 
tential corrugations will in the presence of lateral force or 
shear tend to weaken the attraction when compared to the 
equilibrium case. Neglecting the enhanced friction close 
to the surface, which could lead to chain flattening effects, 
our tentative conclusion is that potential corrugation will 
act similarly as hydrodynamic lift effects and favor poly- 
mer desorption, in line with simulation results for charged 
polymers 21 . Furthermore, since the effective corrugated 
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surface potential becomes more short-ranged compared to 
the equilibrium case, the desorption transition should be 
even more discontinuous compared to a homogeneous sur- 
face potential. 

Interesting effects left out in the present study include 
multi-chain effects, chain stiffness effects and structural 
effects that go beyond modelling chain monomers as a ho- 
mogenous spheres. 
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